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We discuss the 'memory effect' discovered in the 60's by Kovacs in temperature shift experiments 
on glassy polymers, where the volume (or energy) displays a non monotonous time behaviour. This 
^q' effect is generic and is observed on a variety of different glassy systems (including granular materials). 

r—{ , The aim of this paper is to discuss whether some microscopic information can be extracted from a 

quantitative analysis of the 'Kovacs hump'. We study analytically two families of theoretical models: 
domain growth and traps, for which detailed predictions of the shape of the hump can be obtained. 
Qualitatively, the Kovacs effect reflects the heterogeneity of the system: its description requires to 
deal not only with averages but with a full probability distribution (of domain sizes or of relaxation 
times). We end by some suggestions for a quantitative analysis of experimental results. 
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I. INTRODUCTION. THE KOVACS EFFECT 



Systems with slow or glassy dynamics often exhibit non trivial behaviour when temperature changes are applied 
within the glassy phase. Since the system is out of equilibrium, one expects that its properties generically depend 
on the history of the system, an effect that is often called 'memory'. However, this general term embraces rather 
different effects. In the recent spin-glass literature, memory is associated to a two time observable, such as the a.c. 
susceptibility (that depends both on the frequency and on the age of the system) or any other response function. It 
has been shown that after a negative temperature cycle, the a.c. susceptibility recovers the exact value it had before 
the negative temperature jump, hence the name memory. This effect would be trivial if the dynamics was totally 
frozen at low temperature, whereas experiments show very clearly that some noticeable evolution in fact takes place 
■ [1-3]. The same qualitative effect, although not as clear-cut as in spin-glasses, has been observed in many other glassy 



materials (polymers, colloids, ferro-electrics, etc.) [4- 

There is however another well known 'memory effect' that was discovered by Kovacs fourty years ago. This effect 
concerns one time observables, such as the specific volume, or the energy density, etc. and clearly shows that the 
non equilibrium state of the system cannot be fully characterized by the (time dependent) value of thermodynamical 
variables. The procedure followed by Kovacs was the following [9]: first, a reference curve is obtained by quenching 
the sample from a high temperature To to a low temperature T 2 , and measuring the time dependent volume V(t) until 
a time t eq where the system can be considered to be in equilibrium. This defines a volume V eq (T 2 ) = V(t eq ). In a 
"j—; . second step, the sample is quenched again from T to a temperature T\ < T 2 , until a certain time t\. The temperature 
is then quickly raised from T\ to T 2 . The time ti is chosen such that the volume just after the jump reaches the value 
— VeqlTi) - whereas in equilibrium (ti — ► 00) one would have V eq (T\) < V eq (T 2 ). Naively, one expects that 
nothing should happen, since the volume is already at its 'correct' equilibrium value V eq (T 2 ) at the new temperature. 
The volume V(t) in fact shows a non monotonic behaviour for t > ti, first increasing and then relaxing back to the 
q I equilibrium value V eq (T2): 

V(t) = V eq (T 2 ) + AV(t), (1) 



where AV > is the 'Kovacs hump', such that AV(t = tf) = and AV(t — > 00) = 0. Note that the condition 
V(ti) = V eq (T 2 ) (and not V(t^[) = V eq {T 2 )) is chosen as to remove the trivial part of the effect, due to the thermal 
expansion of the fast (local) degrees of freedom. This subtlety in the Kovacs protocol is in fact quite important, as 
will be clear below. 

Fig. 1 reproduces the original results published by Kovacs in 1963 [9], obtained on polyvinyl acetate. The Kovacs 
effect shows unambiguously that other 'internal' variables, beside the volume, are needed to characterize the out of 
equilibrium state of the system, and that these variables did not reach their T 2 equilibrium value at the end of the 
first stage. The memory in this case refers to the fact that these internal variables keep track, to some extent, of 
the system history. To avoid confusion between the different types of memory effects, we will follow [10] and call the 
above phenomenon the Kovacs effect. The Kovacs hump is characterized by its height AVk, and by the time tk for 
which the maximum is reached: AV(t = t\ + tk) = AVk- Qualitatively, the height AVk grows with the temperature 
difference T 2 — T\ (it should obviously be zero for T\ —T 2 ), whereas the time tk decreases when T 2 — T\ increases. 
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FIG. 1. Isothermal evolution at T2 = 30° C of the relative variation of the volume (xlO 3 ) in polyvinyl acetate: after a direct 
quench from T = 40°C* to T 2 = 30°C (1); after quenches from T = 40°C to Ti = 10° C (2), 15°C (3), or 25°C (4) followed by 
rapid re-heating at T2 = 30°C. Data taken from A. J. Kovacs, Adv. Polym. Sci. 3, 394 (1963). 

A similar effect was recently reported in the context of granular materials [11]. In the first stage of another type 
of experiment one 'taps' the system with three different amplitudes -say weak, moderate and strong- during a time 
chosen such as to reach a certain density, identical in the three cases. In the second stage of the experiment, the 
tapping amplitude is chosen to be moderate. The density just after the amplitude 'jump' is recorded. If the state 
of the system was only described by its density, the evolution of the density after the jump should be identical for 
all three situations, and follow the 'moderate' reference curve. This is not the case: as for the polymer glass, the 
weakly tapped system first has to dilate before it is able to resume its compaction, whereas the strongly tapped system 
compacts faster than the reference system just after the jump [11]. 

Finally, the same effect was recently observed in a numerical simulation of three dimensional spin-glasses [12] and 
in a realistic model of molecular liquid [13]. In spin-glasses, the energy density reveals the characteristic Kovacs 
hump when the temperature is raised ; the height of the hump and the time of the maximum behave qualitatively as 
in polymer glasses. Features similar to the Kovacs effect have also been identified experimentally in dipolar glasses 
[14] and spin glasses [15]. Since the Kovacs effect seems to be rather ubiquitous, a natural question is whether 
the underlying physics is the same in all these systems. Stated differently, can the effect select between different 
microscopic models of glassy dynamics? 

The aim of this somewhat didactic paper is to discuss some simple models that allow to shed light on the above 
questions. In these models, the 'internal' variables referred to above appear as a whole distribution function (of 
domain sizes, or of relaxation times) of which only the mean is fixed by the experimental protocol, whereas the shape 
of the distribution keeps track of the system history. We show that the Kovacs effect is indeed rather generic, but that 
the detailed shape of the 'Kovacs hump' could reveal some useful microscopic information on the underlying glassy 
dynamics (see also the discussion in [16]). We first discuss models where slow dynamics is due to a coarsening mech- 
anism, and recall and generalize the main results of [10]. We then turn to the Kovacs effect in the trap model, where 
detailed calculations can be performed. We end the paper with some suggestions for further analyzing experimental 
results, with the hope that the Kovacs effect could help identifying distributions of relaxation times, and/or provide 
some indirect evidence for a growing length scale in glassy systems. 

II. THE KOVACS EFFECT AND DOMAIN GROWTH 

The simplest out of equilibrium system is the one-dimensional Ising model with Glauber dynamics. This system does 
not order at any non zero temperature, but at sufficiently low temperatures the equilibrium domain size £ becomes 
large and for times shorter than the equilibration time, the dynamics is governed by the growth of the typical domain 
size as the square root of time. The energy, which is simply related to the average density of domain walls, plays in 
this model the role of the volume in Kovacs' experiments. When the system is prepared at T\ for a time t\ such that 
the average distance between the walls is equal to the equilibrium size at T 2 > T\ , the out of equilibrium distribution 
of domain sizes at T\ is more sharply peaked around its mean than the corresponding equilibrium distribution at T 2 
-see Fig. 2. In particular, the number of small domains is depleted from its equilibrium value. Upon heating, the 
first effect is that some extra domain walls nucleate within the larger domains, causing the number of small domains 
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(and the energy) to increase. The exact shape of Kovacs' hump can be computed in this model [17], and is found to 
be linear in time for small times, with a slope that increases with the temperature difference T2 —T\, before reaching 
the (exponential) relaxation curve describing a simple quench from high temperatures. Note that the relaxation time 
is finite for all T > in this model; the rate of the final decay only depends on T2, but not on T\. As discussed by 
Brawer [17], this is qualitatively similar to the experimental curves reported in Fig. 1. 
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FIG. 2. Distribution of domain sizes in the one-dimensional Ising model corresponding to two different temperatures, such 
that the average domain size is identical in the two cases (from [17]). The out of equilibrium distribution Ti < T2 (dotted 
line) is more sharply peaked than the corresponding equilibrium distribution at T2 (plain line). Upon heating, small domains, 
initially less numerous, quickly appear within large ones. 

In systems where the equilibrium domain size is infinite, or very large compared to the dynamical length corre- 
sponding to the experimental time scale, 1 the mechanism for the Kovacs effect in coarsening systems is more generally 
the following [10,12]: after a time t\ spent at Ti, the system orders up to a scale l\ = l(t\,T\), leading to an excess 
energy (over the bulk contribution) due to the presence of domain walls with typical scale l\. This excess energy 
density behaves as if ~ d , where is the exponent giving the scaling of the excess energy of a domain with its size 
(for example, 6 = d — 1 for the Ising model, and 9 = d — 2 in the XY-model). When the temperature is increased to 
T 2 , the bulk energy density (within the domains) is suddenly too low compared to the equilibrium value at T 2 . This 
bulk contribution to the energy density therefore increases rapidly by nucleating new domain walls within the large 
preexisting domains of size l\. This picture was actually suggested in [14] to interpret an 'overshoot' effect in dipolar 
glasses which, with hindsight, is the precise counterpart of the Kovacs effect in these materials. 

For larger times, the primary coarsening process resumes and the density of domain walls decreases, leading to 
a decrease of the total energy density. This decay is a priori expected to dominate when the length scale ^(t, T2) 
associated to the dynamical processes initiated by the temperature change becomes of the order of £\, i.e. after a 
time r* such that: 

t( T ;T 2 ) = e 1 . (2) 

However, the time tk at which the maximum of the Kovacs hump occurs turns out to be, in general, much smaller 
than t* (but still much larger than the microscopic time scale To). More precisely, the following picture emerges from 
the exact computation of [10]: 

• In the limit where t(t, T 2 ) and l\ are much larger than the lattice spacing a, the fast initial nucleation processes 
have taken place, and one can expect the energy density E to take the following scaling form: 



1 If the equilibrium correlation length is small, then the Kovacs effect is trivial in the sense that only fast degrees of freedom 
need to reequilibrate. Only if the maximum of the Kovacs hump occurs at times much larger than the atomic time scales 
will one observe a non trivial effect when following the Kovacs protocol, where the contribution of fast degrees of freedom is 
removed by choosing E(tf) = E eq (T2)- (See the discussion in the introduction). 
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AE = E(h + t) - E eq (T 2 ) = AE K T 



£(t,T 2 ) 



£(t,T 2 ) > a 



(3) 



where AEk is the height of the Kovacs' hump, and J-(u) ~ u e d when u — > oo. Using the fact that AE should 
not depend on t\ at large times, one finds AE K ~ £f~ d , which means that the energy scale of the hump is of 
the order of the excess energy stored in the domain walls at Ti . As shown in [10] , the above scaling form indeed 
holds exactly for the 2D XY model in the ordered critical phase, for which O = d — 2. One finds in that case 
Txy(u) = (1 + u 2 )- 1 . 

• In the short time limit £(t 7 T 2 ) ~ a, one expects a nucleation contribution to AE responsible for the Kovacs 
hump. In the case of the critical XY model, where the thermal correlation length £ is infinite, one finds a 
power-law contribution [10]: 



AE w AE K 



1 - 



t(t,T 2 ) 



(-)- 



(4) 



Note that this contribution vanishes for t = 0, since £{t = 0,T 2 ) — a, but cannot be written as a scaling function 
of i/i\. This is at the origin of the difference between tk and t* . The above results only hold if l\ <C £• In 
the other limit where the correlation length £ is small, the above power-law is replaced by a fast exponential 
convergence. In this case tk ~ t , and the Kovacs effect becomes trivial (it would actually disappear if the 
Kovacs protocol was used - see the previous footnote) . 

In this domain growth scenario, one finds the length scale £\ (and therefore r*) to be a decreasing function of 
T 2 —T\. Physically, this is because the bulk energy contribution is lower for smaller Ti; the residual domain wall 
energy density (~ if ~ d ) must then be larger in order to ensure that in the Kovacs protocol, the time t\ is determined 
such that: 

£?(ti,Ti) =E eq (T 1 )+lf- d = E eq (T 2 ). (5) 

For small T 2 — Ti, one thus expects a linear relation AEk C(T 2 — Ti.), where C is the specific heat. 

Therefore, the qualitative dependence of both r* and AEk ~ if ~ d with T 2 — T\ is correctly predicted by this picture. 

If the length £(t,T) grows as a power of time, then from Eq. (3), AE/AEk is found to be a scaling function of 
t/r* in the limit of large times, where the initial (non scaling) contribution due to nucleation vanishes. Due to this 
non-scaling contribution, the scaling function T has a non zero value for small arguments: .F(0 + ) > 0. Therefore, 
in the domain growth scenario (including the equilibrium case discussed by Brawer and recalled above), the Kovacs 
hump does not rescale as a function of t / tk, because the position of the maximum tk is determined by the non 
scaling nucleation contribution. The time should rather be rescalcd by t* determined such that the amplitude of the 
hump has decreased by a factor two (say). By the same token, one expects to see systematic deviations from scaling 
in the regime t <c t* , due to the non scaling contribution of nucleation process. 

We now turn to another soluble model that, interestingly, predicts a variety of shapes for the Kovacs hump, which 
in some regimes are very similar to the ones predicted by the domain growth model. 



III. THE KOVACS EFFECT IN THE TRAP MODEL 



A. Definition of the model 



A simple model exhibiting glassy behaviour is the trap model, which has been extensively studied in the literature 
[18-20], and generalized to describe the rheology of soft glassy materials [21], or the dynamics of contacts in granular 
media [22]. In this model, a particle is trapped in potential wells, and can escape only through thermal activation. 
The depth (energy barrier) of the well is a random variable E > with an exponential a priori distribution p{E) = 
Tg 1 e~ E l Ta . When the particle is in a trap j of energy Ej, it will escape after a time At distributed according to 
Pj(At) = tJ e~ At / Tj , where Tj = r e Ej l T is the mean trapping time of the site j, and then chooses a new trap among 
all the others with a uniform probability. The microscopic time scale To is taken as the time unit in the following. 
The energy scale T g turns out to be also the phase transition temperature. For T > T g , the system equilibrates and 
behaves like a 'liquid', whereas for T <T gi the lowest energy states become the most probable ones and the system 



4 



never stops aging. Of course, this model should not be considered as a realistic microscopic model, but rather as a 
coarse-grained phase-space model -see the discussion in [23,24]. Also, an exponential distribution of energies might 
not be the most appropriate description of a given system. For example, recent simulations of Lennard-Jones systems 
[25] have shown that a Gaussian distribution of barriers is in fact more adequate. As noted in [20], the results of the 
exponential trap model can be extended to that case. 

Due to its simplicity, this model allows one to obtain analytic expressions of many quantities of interest. As for 
coarsening models, we have chosen the energy as the natural observable that plays the role of the volume in Kovacs' 
experiments. 

Let us now present the explicit calculation of the energy as a function of time, with the temperature protocol 
defined in the introduction. However, since fast degrees of freedom are absent in the trap model (there is no 'bottom 
of the wells' dynamics), one does not need to distinguish between t± and tf, as is important both experimentally and 
in models with microscopic degrees of freedom (see above for a discussion of this point). Two different cases have 
been considered in details. In the first one, the temperatures T\ and Ti are both above T g , but close to it, so that the 
system eventually equilibrates, but with very long relaxation times. In the second case, both temperatures are below 
T g , so that the system is in the aging regime where equilibration is never achieved. Finally, we only briefly discuss 
the 'mixed' case where T± <T g < T 2 . The original Kovacs experiment corresponds to the first case, since the volume 
is seen to relax towards its equilibrium value at T 2l used as the reference energy. In the second case, the time t\ at 
which temperature is shifted is in fact arbitrary, but interesting scaling properties appear. 



B. Case T > T g : relaxation towards equilibrium 

We shall use a continuous energy description (see [20]), i.e. the system is described by the probability P^(E,t) to 
be in a state with an energy (barrier) E at time t and temperature T, which evolves according to the following Master 
equation: 

^(E, t) = - e- E ^ T P T (E, t) + u{t)p(E) (6) 

with u(t) = J °° dE' e~ E / T P^{E' , t) is the average hopping rate. For T > T g , Pt(E, t) relaxes towards the equilibrium 
distribution P^ q (E) = Z~ x e E / T . So the interesting quantity to study is the deviation from equilibrium, i.e. the 
distribution pr(E, t) defined as pt(E, t) — Pr(E, t) — P^?(E). Let us first focus on a simple isothermal quench from a 
given initial condition Pq(E) = P^? o (E). The evolution of Pt{E, t) can be computed using a time Laplace transform, 
and if T > T, the asymptotic behaviour of the distribution becomes independent of the initial condition Po(E), 
yielding: 

m s) - {P9 \ P l e ^E P)E m r(2 - e) en (7) 

where — 1/T and 6 = T/T g is the reduced temperature. Let us define the energy deviation £t(*) = \Er(t) ~ E^\. 
Note that in the following, energies are understood to be true physical energies, i.e. the opposite of the energy barriers: 
£ T(t) = — Jo°° dE E pt(E, t). This last quantity can be computed from p(E, s), which gives: 

M*) = ^[T(0)lnt-r'(0)]. (8) 

Hence, the energy relaxation above T g is (up to a logarithmic correction) a power law with an exponent that becomes 
small for T — > T g . The time t\ when the temperature has to be raised from T\ to T2 in the Kovacs procedure is 
defined by E^ x (t\) — E^, or equivalently £Ti(^o) = — E^ . Thus t± is determined by the equation: 

-J^^Olnii-m)],,^^. (9) 

Note that in order to be consistent, the above equation assumes that 9 2 — 61 -C (6*i — l) 2 <C 1, in which case 
h >r (= 1). 

Now using the distribution px 2 {E, to) = PT X (E, t ) + P^? (E) — P^ (E) as initial condition in the Master equation, 
one can compute the further evolution of the energy at Ti at time t\+t. A time scale r* = t[ naturally appears (with 
7 = 61/62). Defining the energy variation AE(t) = E(t\ +t) — E(t\), one finds in the short-time regime 1 <C t <C r*: 



5 



AE(t) 



Ti 



1 



1 \ T 2 T 2 

1 + pr^r ^177 t ln< + ^ + ^r[ lnt + Tel 



(10) 



where — — r'(l) is the Euler constant. Interestingly, this behaviour is very similar to that found for the coarsening 
process. Indeed, in the limit 9 2 — 6\ <C {0\ — l) 2 , the two power laws in the previous equation, 9 2 — 1/7 and # 2 — 1 
are very close to each other, and the expression can be simplified as: 



Therefore, the maximum of the Kovacs hump is given (in the considered limit) by: 

AE K 



t 



(0i 1) 



2 ' 



(11) 



(12) 



The approach to this maximum is described by a power law of time with a logarithmic correction. This is not very 
different from the coarsening model discussed in the previous section. Note that the height of the hump is again 
linear in T 2 — T\ for small temperature differences, as was the case for domain growth. Note also that r* = t[ is a 
decreasing function of T 2 — Ti , in agreement with experimental results. 

In the long time regime f>T*, one recovers as expected the isothermal quench result Eq. (8) at temperature T 2 : 



T 2 



AE(t) = lni-r'(0 2 )] 
This late time result can again be put in a scaling form (up to logarithmic corrections): 

AE(t) = AE K g(j^j r*=tj 



(13) 



(14) 



but the early time regime Eq. (11) fails to scale. Only when Ti,T 2 — > T g , does one find that the maximum time tk 
coincides with r*. More generally, and as for domain growth, one has tk <Cr*. 



C. The aging case (Ti, T 2 <T g ) 

We now turn to the aging case where the shape of the hump is found to be qualitatively different. We consider the 
case where both T\ and T 2 are less than T g . In this case, the system never converges to an equilibrium state, but keeps 
on aging, so that the situation is different from that of the Kovacs original experiment, but could in principle also be 
investigated experimentally. Since the equilibrium energy at T 2 does not exist, we choose to shift the temperature 
from Ti (initially reached at t = 0) to T 2 after a waiting time t w (which plays the role of t\ in the previous sections). 



1. Probability distribution and Green function 



The continuous energy Master equation (6) does not admit anymore a stationary solution. The resulting dynamical 
distribution can be computed using the Laplace transform Pj(E,s w ) (in the time domain) of Pj(E,t w ), where t w is 
the waiting time since the quench from high temperatures. One finds, in the asymptotic regime s w — > (or t w — > 00): 

Pr(E, s w ) * Ut(E, S )^^ {1 + a ^E ){aeP E r ( 15 ) 

where 9 = T/T g . Since stl^(E,s) is a function of the product se l3E , Il^(E,t) depends only on the scaling variable 
£ = el iE jt. Then one can turn to the computation of the energy variation after the temperature shift. 
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FIG. 3. Plot of AE(t w ,t) in the trap model as a function of the scaling variable t/tZ for t w = 10 5 , 10 6 , 10 7 and 10 s (Monte 
Carlo data), 6 1 — 0.5 and 9 2 = 0.6. Inset: comparison between Monte Carlo data (t w = 10 8 J and the analytical prediction of 
the short time behaviour -see Eq. (21 )- for the same temperatures as above. 



2. Computation of the energy variation 

The detailed calculation of the energy variation AE(t w ,t) between time t w (when temperature is shifted from T\ to 
T 2 ) and t w +t is given in Appendix A. Here we shall only summarize the main steps of the calculation, and emphasize 
physical interpretations and conclusions. From a technical point of view, a useful tool in order to compute AE(t w ,t) 
is the Green function Gt(E, Eq, t) defined as the probability to have the energy E at time t w +t given that the energy 
was E at time t w . This Green function is computed, as for Pt(E, t w ), using a Laplace transform, with s the Laplace 
variable. One finds, for sr <1 (t> r ), the following asymptotic expression: 

P P E o 1 
Gt(E, E ,s) = tt -^ S(E Eo) + -^—— m - U T (E, s) (16) 

Thanks to the Markovian properties of the dynamics, the Green function does not depend on t w , but only on the 
time difference t. One can express the average energy E(t w + t) at time t w + t using the Green function: 



E(t w + 1) 



poo POO 

= - dE dE EG Tl (E,E ,t)P T2 {E ,t w ) (17) 
Jo Jo 

where the minus sign accounts for the fact that the variable E (i.e. the energy barrier) is actually the opposite of the 
true energy. The energy variation AE(t w ,t) is then: 

AE(t w ,t) = E(t w +t)-E(t w ) (18) 

pOO pOO 

= -/ dE dE (E-E )G Tl (E,E ,t)P T2 (E ,t w ) (19) 
Jo Jo 

After a few calculations (see Appendix A), one can show that AE(t w ,t) exhibits a kind of 'sub-aging' scaling (see 
also [26]): 

AE(t w ,t)=il>(^j, (20) 

where 7 = T\jTi < 1. One sees that the energy evolves on a typical time scale given by tk — t* = t^, which is 
expected from a simple activation argument. One can also study the asymptotic (short time and late time) behaviour 
of this scaling function. Concerning the short time behaviour, it is necessary to distinguish between two cases. 



If 7 > 1 — 6*i (small temperature shifts) then AE(t w ,t) is found to be singular at short times: 

/ t \ (1-0O/7 

A£(^,t)~^> ^ t«C (21) 
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• If on the contrary 7 < 1 — 61, one finds a linear t dependence in the short time regime (with logarithmic 
corrections) : 

AE(t w , t) ~ K< fin *f + t « tl (22) 



The coefficients K > , K< and C appearing in Eqs. (21, 22) are given in Appendix A -see Eqs. (64, 70, 71)- and 
are found to be positive for Q\ < 2 . Therefore AE(t w ,t) is positive for short times, and the Kovacs effect has the 
expected sign. 
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FIG. 4. Plot of AE(t w ,t) in the trap model as a function of the scaling variable t/tZ for t w = 10 6 (Monte Carlo data), in the 
case of a negative temperature shift: 6\ = 0.6 and 62 = 0.5. Inset: comparison between Monte Carlo data and the analytical 
prediction of the short time behaviour -see Eq. (21); note that finite size effects are strong in this case. 



Note also that the coefficient K > vanishes linearly when 6\ — > #2- This is expected: if no temperature jump occurs, 
the energy variation should be regular, i.e. linear in t. Moreover, if Q\ > 8 2 (negative temperature shift, 7 > 1), the 
above calculation is still valid, with a negative K > , and a non trivial exponent (1 — 81) /j. The Kovacs hump becomes 
in this case a Kovacs trough. Monte-Carlo data are compared with these analytical predictions in Figs. 3 and 4, 
showing a rather good agreement. Note that the scaling Eq. (20) is only approximate for finite t w . A better rescaling 
can be obtained in the case 0\ < 62 by plotting AE(t, t w )/AEx as a function of t/rK , where AEk is the maximum 
value of AE(t w ,t), reached at t — tk — t^. Eq. (20) means that asymptotically, AEk becomes independent of t w . 

Finally, the long time behaviour is easy to analyze: one can show that Pt(E, t w + 1) behaves asymptotically in the 
same way whatever the initial condition P^(E 1 t w ). The system behaves, at late time, as if it had been quenched 
directly from high temperature (see Fig. 5). This means that in this limit E(t w + 1) docs not depend on t w , but only 
on t: 

E(t w +t)~ E late (t)=T 2 [V '(l)-7r cot n6 2 }-T 2 lnt t » t£ (23) 

where Ei ate (t) is the average energy at a (large) time t after a quench from a high temperature. So AE(t w ,t) is 
simply given by: 

AE(t w ,t) = E late (t) - E(t w ) t » tl (24) 
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FIG. 5. Comparison of the late time behaviour of E(t w + t) = AE(t w , t) + E{t w ) in the trap model as a function of t, for 
temperatures Oi = 0.3 and 0.4, O2 = 0.5 and waiting times t w = 10 7 and 10 6 respectively (dot-dash) with a direct quench from 
infinite temperature to 82 = 0.5 (t w = 0, full line). The asymptotic analytical prediction is also shown (dashed line). 



D. The mixed case (Ti <T a <T 2 ) 

For completeness, and since this is also an interesting situation, we briefly mention the results obtained in the trap 
model in the case where the glass transition temperature T g lies between T\ and T 2 : T\ < T g < T 2 . This case is 
worth studying, since the system can eventually equilibrate at the final temperature T 2 , with long relaxation times 
(assuming T2 is close to T g ), but T\ can be varied in the whole range < Ti < T g and not only in the vicinity of T g . 
Interestingly, one finds the same short time singularities as in the aging case (Ti <T 2 <T g ) studied above : 

/ t \ (l-0i)/7 

AE(t u t)~^J 1-0 1<7 (25) 

AE(t u t)^^ 1-0 1>7 (26) 

with however prefactors and logarithmic corrections which are different from that found in the aging case. In the long 
time regime, one naturally finds a convergence of AE(t,t\) proportionnal to i - ^ -1 ), as in the case T g < T\ < T 2 . 
As a result, one sees as could have been expected that the short time regime is generically dominated by the thermal 
history before the temperature shift, whereas the long time behaviour depends only on the final temperature T 2 . 



E. Discussion — 'Fronts' in the energy distribution 

It is interesting to discuss how the distribution of energies P(E, t) evolves when the temperature is shifted from T\ 
to T 2 > T\. This is illustrated in Fig. 6: at the lowest temperature, the probability of small (negative) energies is 
depressed. When the temperature is raised, the system obviously re-equilibrates fastest in the region of small energies 
since this corresponds to the smallest relaxation times. The probability 'hole' is thus rapidly filled, leading to an 
increase of the average energy, and thus to the Kovacs effect. As time increases, the equilibration progresses as a 
kind of 'front' in energy space, as shown in Fig. 6. Only at later times does the peak of the distribution move to 
larger (negative) energies. It is interesting to realize that this picture is in fact very close the one emerging from the 
coarsening model where short scales re-equilibrate fast and lead to an increase of the average energy, before larger 
length scales resume the coarsening process (see Fig. 2, and the discussion of Section II). 

The conclusion from the 'domain growth' interpretation of the Kovacs effect presented in the previous section 
is that a quantitative analysis of the Kovacs effect might give one a unique tool to investigate experimentally the 
problem of growing length scales or the statistics of trapping times in glassy systems, a topic of huge current interest 
[27-30,25,31,32]. However, as demonstrated above, the quantitative predictions of the trap model are in fact very 
similar to that of domain growth. As discussed recently in [33], the physical difference between the two pictures is 
not as obvious as it might first seem. In particular, the trap model description implicitly assumes the existence of an 
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underlying 'coherence length' [23]; conversely, domain growth models may naturally generate a non trivial distribution 
of relaxation times [33]. A possible discrimination might lie in the temperature dependence of the short-time and 
long-time exponents that describe the Kovacs hump. While a temperature dependence is expected in an activated trap 
like description, it is less natural for power-law domain growth. On the other hand, more complicated (logarithmic) 
growth laws can mimic power-laws with a temperature dependent exponent [34,12]. 




FIG. 6. Dynamical energy distribution P(E,t w ,t w +t) at time t w +t after a (positive) temperature shift at t w , witht w = 10 4 , 
0i = 0.5 and 62 = 0.8 (trap model). Time t ranges from t = to t = t w . One clearly sees the propagation of a 'front' at small 
energies, associated to the re-equilibration of the short time scales, before the global drift of the distribution starts again. 



IV. SUMMARY AND CONCLUSIONS 



Although the Kovacs effect has been known for fourty years, its quantitative interpretation has not been much 
developed until recently. In view of the fact that this effect is generic and observed in a variety of different 'glassy' 
systems and models (such as the ones studied in the present paper), it is important to establish which type of 
microscopic information one can extract from the quantitative analysis of the 'Kovacs hump'. Qualitatively, the 
Kovacs effect reflects the heterogeneity of the system: fixing the overall (macroscopic) value of the volume or energy 
does not prevent the existence of local fluctuations, which keep track of the system history. A more complete 
description of the system therefore requires to deal with full distribution functions, and not only with averages. In 
the two models studied in this paper, this distribution function is that of domain sizes in domain growth models, and 
that of relaxation times in the trap model. These models lead to precise, quantitative predictions for the shape of the 
Kovacs hump, which are summarized in Table I, and by the following phenomcnological equation that describes the 
Kovacs hump [inspired from both domain growth and trap models -see Eqs. (3,11,13)]: 



AE = AE K 



(t + t*) v ¥ 



(27) 



In domain growth models, the exponent v is equal to (d — @)/z, where z is the dynamical exponent relating length 
and time: I ~ t x l z . The first term corresponds to the long time contribution of already grown domains, and the 
second term to the excess energy created by the transient nucleation of new domains. A similar interpretation of the 
two terms can be given within the trap model, with v = T/T g — 1 (for an exponential distribution of barriers). Note 
that from Eq. (27), one sees that the data should re-scale as AEk £/(-tt) when r* becomes large. The position of the 
maximum t k extracted from Eq. (27) is, in the limit t* » 1, given by: 

t k « (w*)^ < t* (28) 

so that the correct scaling variable is not t/rK, except in the limit v — > 0, where tk and r* coincide. 

We hope that these results will motivate new, systematic experiments. It would in particular be valuable to test 
Eq. (27). It would be very interesting to extract from a detailed analysis of the Kovacs effect a quantitative determina- 
tion of the distribution of relaxation times in glassy systems, and its temperature dependence. This distribution could 
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then be compared with other direct, dynamical determinations. Another situation worth investigating experimentally, 
suggested by the present study, is the out of equilibrium (aging) Kovacs effect, where both temperatures are kept well 
below the glass temperature T g . 

Finally, from a theoretical point of view, it would be worth studying the predictions of the mean field (p-spin) 
spin-glass for the shape of the Kovacs hump. As is well known (see e.g. [35]), the dynamical equations for this system 
are identical to the Mode-Coupling equations for structural glasses. Although we expect on general grounds that the 
results of this model should be again quite similar to those obtained in the present study, it would be interesting to 
check this assertion in more details. 





Domain growth (£ eq — oo) 


Trap model Ti > T g 


Trap model Ti < T g 


Preparation time ti 


^i,Ti) w - d ocT 2 -Ti 


t 1 -" 1 lnii cx T 2 - Ti 




Height of the hump AEk 


T 2 - Ti 


T2-T1 


T2-T1 


Characteristic time r* 


£(t*,T 2 ) =£(ti,Ti) 


r*=t 1 1 l/ ' 1 ' 2 


T*-t ri/ ' r2 

l - o w 


Hump time tk 


TK ~ (r*)^TT < T * 


T K ~ (r*)^TT < T * 




AE at early time 


AE K (1 - t°- d ) 


AE K (1 -t-"\nt) 


AE K (t/T*) { - 1 - til)/ ~< 


AE at late time 


AE K (li/£) a - a 


AEk{t* /t) v 


-T 2 \n{t/tZ) 


Exponent v 


u = {d-Q)/z 


V = 02~l 





TABLE I. Summary of the different results and regimes for the Kovacs hump, in the limit where Ti — + T 2 7. We denote by 
ti or t w the time spent at the lowest temperature Ti, and z the energy and dynamical exponents for domain growth, 7 the 
ratio 7 = T1/T2 and Q\ t2 the reduced temperatures #1,2 = Ti j2 /T g . 
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APPENDIX: DETAILED CALCULATION IN THE TRAP MODEL 



A. Probability distribution and Green function 

The Master equation of the trap model reads: 

'-(E,t) = - e - 0E P T (E)+w(t)p(E) 



dP T 



dt 



(29) 



with uj(t) = / °° dE'e l3E P^(E' ,t), and (3 = ^. E is a positive variable, the depth of the traps, and is actually the 
opposite of the true energy of the states. This Master equation has to be supplemented by an initial condition: 



P T (E,t = 0) = P (E) 



(30) 



where Pq{E) is a given (arbitrary) probability distribution. We also take an exponential density of states, p(E) 
Tg 1 e~ E / Ts . Introducing the Laplace transform Pt(E, s) with respect to t defined as: 



the Master equation becomes: 



Solving for Pt(E, s), one has: 



Pt(E,s) = / dte- st P T (E,t) 
Jo 



sPt(E, s) - P (E) = -e- fiE Pr(E, s) + <2>{s)p{E) 



Pt(E,s) 



P (E) + £>(a)p(E) 



s + e~f 3E s + e~P E 



(31) 



(32) 



(33) 



u)(s) is determined by multiplying Eq. (33) by e & E and integrating over E. The distribution P<r(E, s) is then given 
by: 



e 0E 1 e l3E o(E) T r 



dE 



P{E) 



with (p(s) defined as: 



f 

Jo 



dE 



ME) 
1 + se^ E 



(34) 



(35) 



Integrating Eq. (34) over E allows to check that Pt{E, s) is well normalized, i.e. J °° dE Pt(E, s) = 1/s. In order to 
compute the variation of the energy after a temperature shift, one has to introduce the Green function Gt(E, Eo,t) 
defined as the probability for the system to have energy E at time t w + 1 given that the energy was Eq at time t w , if 
the bath temperature is T. Note that since the process is Markovian, the Green function depends only on the time 
difference t, and not on t w . The Green function in Laplace space Gt(E, E , s) is straightforwardly obtained from 
Eq. (34) choosing P (E) = 5(E - E ): 



G T (E,E ,s) - 



l + se 0E o S ( E E 0)+ s l + se (3E l + s [js 



O 0E 



p(E) 



dE 

l + seP E 



(36) 



As shown in [20], the energy distribution P^{E 1 t) takes a scaling form for large times. Indeed, from Eq. (34), one 
has for s — > 0: 
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P T (E,s) 



1 e^ E p(E) 
s 1 + s eP E 



dE 



e» E p{E) 



(3e 



0E 



7T (l+se^)(se^) s 



U T (E,s) 



(37) 



which defines the asymptotic distribution U.t(E, s). The reduced temperature 9 = T/T g has also been introduced. 

E /T 

Its inverse Laplace transform II x (E, t) satisfies a scaling relation in the variable £ = : 



One finds g(£) by inverting the Laplace transform given by Eq. (37): 



sin7r# 1 
^ ) = ^)? e " 

So for large times, the Green function is given by: 

P I3E„ 



/■V€ 
^0 



G t(^, ^o, «) = 1 + se , Eo S(E E ) + 1 + sefjE(> U T (E, s) 



(38) 



(39) 



(40) 



We now consider the following thermal history: at the initial time, the system is quenched from To > T g to T\ < T g \ 
at time t w , it is re-heated to a temperature T 2 satisfying T\ < T 2 < T g . One is interested in the subsequent evolution 
of the energy, at time t w + 1. The probability to have energy E at time t w + t, given this thermal history, is: 



/>oo 

P(E,t w ,t) = / dE G T2 (E,E ,t)P Tl (E Q ,t w ) 
Jo 

Taking the double Laplace transform with respect to t w and t: 

t>OC 

P(E,s w ,s)= / dE G T2 (E,Eo,s)P Tl (E ,s w ) 
Jo 

Using the asymptotic expressions Eqs. (37) and (40), one can write: 



(41) 



(42) 



P(E, s w ,s) = TT -^ n Tl (E, s w ) + n T2 (E, s) r dE n ^ f,E 



(43) 



B. Evolution of the average energy and scaling relation in the aging regime 

Bearing in mind that the energy of a given state is the opposite of the energy barrier E, the mean energy E(t) at 
time t is defined by: 



-E(t) = r 

Jo 



dEEP T {E,t) (44) 



Taking into account the thermal history introduced in the preceding section, we define the energy variation between 
time t w and t w + t: 

AE(t w ,t) = E(t w +t)-E(t w ) (45) 
Computing the double Laplace transform yields: 

/>OC />oo 

-AE(s w ,s)= dE dEo(E-E Q )G T2 (E,E ,s)P Tl (E ,s w ) (46) 
Jo Jo 

Expanding this equation, one finds for s w to and sto <C 1: 

- AE(s w ,s) = J~*EE n T2 (E, s) [ dE, - I J™ dE E (47) 

= i{s)J(s w ,s)--k(s w ,s) (48) 
s 
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where I(s), J(s w ,s) and K(s w ,s) denote respectively the three integrals appearing in Eq. (47). Making the change 
of variable r = e@ 2E in / and r = eJ 3 ^ in J and K, one has: 



1(a) 
K(s w ,s) 



sin7T# 2 
sin 7r^i 



dr- 



T 2 lnr 



(1 + st)(st) 6 
dr 



7T (1 + STT)(1 + Sw r)( 



sin tt9i 



dr Ti In r 



7T (l + ST-r)(l + s M ,T)(s tu r) e i 

with 7 = Ti/T 2 ; I(s) can be computed using the identity lnr = dT a /da\ a=0 ] one finds: 

T 2 

I(s) = (7TCOt7T02 — Ins) 

Let us show that AE(s w ,s) satisfies a scaling relation. Note first that for s — ► 0, J(s w , s) is of the form: 

/(s^/s 1 ^) 



— f 

^ Jo 



du ■ 



l + u"t 



where f{x) = (sin7r0i)/[7rx 6 ' 1 (1 + x)], and u = s x ^t. In the same way, K(s w ,s) reads: 



K(s w ,s) = -T 2 lnsJ(s w ,s) + J du^-^ 



In i 



(49) 
(50) 
(51) 

(52) 

(53) 

(54) 



Coming back to AE(s w , s), one has from Eq. (48) that terms in Ins cancel, and one gets: 



-AE(s w ,s) 



1 



s l+l/7 



7rT 2 COt 7T# 2 



Jo l + « 7 Jo 1 



/(s^/s 1 ^) 



ln?i 



,i*h (55) 



which implies a simple scaling form AE(t w ,t) — ^(t/t^). This is easily shown by computing the Laplace transform 
of this scaling form: 



(56) 



Let us make the following changes of variable: t — xf^ (at fixed t w ), and then t w — v/(sx) 1 ^' 1 (at fixed x). One finally 
gets: 



1 



,1+1/7 



™V>(*)/ d^exp(-^^-^)= * ^(-^) 



(57) 



which indeed gives back the expected scaling form in Laplace space. 



C. Short time behaviour 

In this section, we shall focus on the short time behaviour of AE(t w ,t), characterized by t <C Q,, or equivalently 

s 3> sj. Note however that we consider only times that are large compared to the microscopic time scale: t,t w ^> 
ro = 1 (s, s ffi < 1). From Eq. (55), one sees that two integrals have to be computed: 

sin7r#i f°° du 

{) = ^^J (l+u->)(l + Xu)(Xu)^ (58) 
sinTr^i [°° \nudu 

B{X) -^^i (1 + U 7)(l + A«)(A«)^ (M) 

where A stands for the ratio s^/s 1 / 7 . In the case A<1, these integrals reduce to: 
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sin7r#i Z" 00 du sin7r^i f°° \nudu 

A{x) = ^x^J JTT^W ( ) = ^a^7o (i + ^K 1 (60) 

on condition that 6\ + 7 > 1. The opposite case, 81 + 7 < 1, will be considered later on. The integrals A(X) and 
B(X) are readily calculated using the following identities: 

f°° dv 7T f°° \livdv 7T 2 C0S7T/Z 

J vv(l + v) sin 71^ 7 1^(1 + u) sin 2 71^ 
Altogether, one finds for AE(s w , s): 

a / \ T27rsin7r0i r 7r „ . „ . 1 / s 1 / 7 ^ 1 . . 

-A£'(s w , .) * 78in ^ (1 _ gi) [cot -(1 - ft) + cot ^]-^ ( — J (62) 

The short time behaviour of AE(t w ,t) is obtained by inverse Laplace transform, in the case 9\ + 7 > 1: 

/ t \ (i-eo/7 

A£(i w ,i)~#> t«t 7 (63) 



where the coefficient if> is given by: 



T 2 7rsin7r6»i[cot^(l - ft) + cot tt6 2 } 



K> = : — rr fl — 1 (64) 

7 sm[^(l-0 1 )]r(0 1 )r(i±^) 

Note that in spite of the minus sign in the r.h.s. of Eq. (63), AE(t w ,t) is indeed positive at short times for ft > #1, 
showing that the energy has to increase first before reaching lower values. However, this coefficient vanishes for 62 = ft 
(i.e. no singularity occurs if temperature is kept constant), and becomes negative for ft < ft. 

In the opposite case, ft + 7 < 1, another approximation has to be used. Making the change of variable v = Xu in 
A(X) and 5(A) -see Eqs. (58,59)- one finds: 



sin7rft f°° dv 
{) = ^X~J [l + (v/\)-r]{l + v)v°i 



= sm^rft /« lnv~lnX 
V ' ttA J [l + (v/X)~t](l + v)v i v ' 

In the small A limit, (v/X) 1 » 1, so that A(X) and B(X) reduce to: 

. ... sin tt9i f°° dv sin7rft f 00 , In v — In A ,„„. 
A(X) ~ x1 / , fl T B(X) ~ -: / dv—jr— (67) 

which are indeed convergent since 7 + #1 < 1. One then finds for AE(s w , s): 

sin7r#i r s i 1 /s 1 / 7 \ 1_7 

-AE(s w ,s) * — cot .ft - ^ cot .( 7 + ft) + Ti In (— J (68) 



The inverse Laplace transform yields: 



t \ t 



AE(t w ,t) ~ K< ( C - In 1 -y i«C (69) 



where and C are given by: 



K< = =7 ^ Sin7r ^ ^ (70) 



r(l - 7 )sin7r( 7 + 6li) 

r'ci - 'vi 

C = 7^ y -r'(2) + 7rcot7r6i2 - tt cot 7r( 7 + ft ) 

L(l - 7) 



(71) 
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D. Long time behaviour 



One can also study the long time behaviour t ^> t%, , which happens to be easier to handle than the short time one. 
Coming back to the starting equation (47), the limit s<sj simplifies a lot the equation, and one gets: 

rp poo pOO 

-AE(s w ,s) = —[ircotir02-ln S } dE U Tl (Eo,s w ) — / dE E fl Tl {E , s w ) (72) 
s Jo s Jo 

S v ' 

l/s w 

The second term is nothing but Cu w E(t w ), which also appears in the left hand side of the equation, due to the 
definition of AE(s w , s). In other words, the long time behaviour of E(t w +t) appears to be the same as if the system 
had been quenched from high temperature to T 2 at time t w . One finally finds: 

AE(t w ,t) = E{t w +t)-E(t w ) (73) 
— T 

E(t w +t) = -Ct 1 — (7TC0t7r<9 2 - Ins) (74) 
s 

= T 2 [r'(l)-7rcot7r6»2]-T 2 lnt (75) 

showing that E(t w + t) is indeed independent from t w and from T\ for times t ^> V w . This result can also be found 
directly without using this particular thermal procedure: if one computes the probability distribution P(E, t) for large 
times t, starting from an arbitrary initial distribution Pq(E), it appears that the asymptotic (large t) distribution 
does not depend on Pq(E): 

E(t w +t,t w )=E late (t) (76) 
where Ei ate (t) is the average energy at a large time t after a quench from high temperature. 
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